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We examine the properties of a quantum reflection trap when particle-interation is included. 
We explore the influence of the particle-interaction on the trapping for different regimes: repulsive 
particle-interaction and attractive particle- interactions in its stable and unstable limit. With varia- 
tional techniques, we calculate the phase-diagram of the quantum reflection trap and determine the 
stable and unstable regimes of the system. 
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^ . Recent progress in experimentally controling BECs [J, [3| has shown, that not only single atoms, but condensed 
^_ ' atomic clouds, Na in these cases, can be quantum reflected as a whole by atom-surface potentials when the kinetic 



energy of the incident cloud lies within the threshold region of the potential. Armed with these results, we suggested 
to trap ultra-cold atoms solely by quantum reflection in [3] . The advantages of quantum reflection over other trapping 
mechanisms are obvious: stability of the reflection behaviour in the threshold region, no need for external auxiliary 
fields to form the trapping potential, because atom-surface interaction is delivered by nature free house. In Q we 
have shown, that the surviving particle density inside the trap as a function of time gives reasonable good results to 
be promising for future investigation. We have shown, that due to atom-surface interaction, an enhancement of the 
' surviving particle density up to 50% is achievable. The trapping times lie somewhat around 0.5 s, excluding mass 
factors of the atomic species used. 
^ ' The model we proposed in Q relies solely on the threshold properties of an atom-surface potential. The trapping 
O^. mechanism, quantum reflection, exploits the fact that cold atoms are reflected by an attractive potential tail without 
' reachinga classical turning point. The mechanism of quantum reflection has enjoyed a growing interest, experimentally 
^ ; 0, i, S 0, 1 and theoretically 0, [13, IHI ■ 

^ ■ The equivalence between the true atom-surface potential and the step-potential at threshold comes from the fact, 
■<!;;j- ' that both obey the same law of reflectivity, see, e.g., Very recently, Madronero and Friedrich, [13], have shown 
T-H , that the dynamics of a wave-packet governed by a step-potential is indeed qualitatively similar to the dynamics of a 
' wave-packet governed by a power-law atom-surface potential. Quantitatively, the differences are negligibly small and, 
as a matter of fact, the step-potential model even underestimates the results of the power-law potential. For quantum 
T-H I reflection, the use of a step-potential as a model for the atom-surface power-law potential has thus been justified once 
more. 

25 ■ Here, we will examine the properties of a spherical quantum reflection trap, which is the simplest choice of a 



trapping system when particle interaction is included. To model our potential in this case, we make the following 
assumptions on the behaviour of the atom-surface potential inside a sphere: 



.1-.^W = -^F^' , (1) 

where L is the radius of the sphere. The assumptions Eq. Ilj are certainly true if the radius of the sphere is, by orders 
of magnitude, larger than the extension of the atomic wave-packet inside the sphere. If this condition is fuUfiUed, the 
inside wall of the sphere, for an atom close to it, must locally look like a plane, In good agreement with these 
requirements, the atom-surface potential inside a sphere can be modeled as a step-potential 

U{r) = - ^ b-^e [r~L], b = (3^ . (2) 



II. THE MODEL OF THE QUANTUM REFLECTION TRAP 



The system of interest is given by the Gross-Pitaevskii equation 

ih^^^ir,t) ^ ^L\/H{r,t) + ^Uir)^ir,t) + 8TTa,^,\^{r,t)f ^{r,t)\ , (3) 
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where the particle density is normaHzed to unity. The scale of the system is naturally given by the spatial extension 
of the trap, the radius L of the sphere. The scaling variables are thus given by 0], 

L p4 2 m iy^ 

The last term in Eq. ^ restores normalization to unity. The atom-surface potential Eqs. lU [5]) shows only a 
dependence in the radial direction, so the angles can be separated off. Assuming the system to be in an s-wave state, 
the scaled wave-function is 

X V47r 
With Eq. ([5]), the angular parts can easily be integrated out which yields for the interaction energy 
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The radial coupling-constant is then given by 

8 TT aint 2 aint 

The radial Gross-Pitaevskii equation describing our quantum reflection trap in scaled form thus finally reads 

i-^^{x,T) = -^i^{x,t) ~ a^^[x-l]^(a:,r) +7 !^^^; ^1 V(x,r) . (8) 

To simplify our notation, we have dropped the subscript of the radial wave-function and the interaction-strength 7. 
The step function on the right hand side in Eq. ([5]) is again the model of the atom-surface interaction. 

In table (|T| we have listed some values of the radial coupling constant for the alkali atoms typically used in BEG 
experiments. Table ([J) shows that the natural influence of the coupling due to the particle-interaction is rather small. 
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TABLE I: Comparison of the scaled potential strength cr and the coupling-constant Tfor some alkali atoms. The radial extension 
of the trap is allover L = 4.47 x 10^ [a.u.] . The pot ential data was taken from [10|, the data for the scattering lengths were 
taken from ^ for Li, from [IE] for Na, and form ^ for Rb . 

such that the results published in Q, where trapping was considered neglecting particle-interaction, are confirmed 
to apply for cases where the particle-interaction strength is not tuned. On the other hand, a tuning of the coupling 
constant by the help of a Feshbach-resonance delivers a wider range for 7, which we now will examine theoretically. 



III. WAVE-PACKET DYNAMICS WITH PARTICLE INTERACTION 



To explore the dynamical properties of an atomic wave-packet with particle-interaction inside a quantum reflection 
trap, we solve Eq. ([H]) with the initial condition 

^(a;,r = 0) = TVx exp[-aa;]6'[l - x] , (9) 

where N is the normalization constant and a is the diffuseness of the wave-packet. The wave-packet Eq. Q provides 
the simplest possible choice for an initial state supporting the main parts of particle density around x — 0, j3|. 
To be more explicit, we take the example of Na-atoms, which have been used in the recent quantum reflection 
experiments with BECs [H, 0. A diffuseness of a = 5, together with the atomic mass of sodium and a trap radius 
of L = 4.5 X 10^ [a.u.] gives an initial kinetic energy i?kin — 1-5 x 10"^^ [a.u.] that corresponds to temperatures of 
approximatly 1 nK. 
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The observables of our examinations are the surviving particle density ps{t) inside the trap, and the scaled energy 
£(t) of the system, as functions of time: 

Ps{r) = / dx \ip{x,T)\'^ , 8{t) = i / dxip*{x,T)dripix,T) . (10) 
Jo Jo 

Numerically, we solve our differential equation ([8|) by using the Crank-Nicholson algorithm and employ absorbing 
boundary conditions, see [l^, to simulate outgoing waves, that have left the spatial region of the quantum reflection 
trap. 

A. Dynamics with repulsive particle-interaction 

In cases where the coupling constant 7 is larger than zero, the particle-interaction is repulsive. The additional 
positive energy due to the self-interaction increases the total energy of the system. An increase of the systems's total 
energy slightly accelerates the decay of the surviving particle fraction inside the trap, because even low-energetic 
components of the wave-packet gain additional energy, which facilitates the escape from the trap. As an example we 
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FIG. 1: Surviving partile densities for ct = 20 and values of 7 = 0,0.1,0.5,1.0,5.0 from top to bottom. The presence of the repulsive 
particle interaction reduces the trapping effect of the quantum reflection by its additional contribution to the total energy of the system. 

have shown in Fig. ([1]) the surviving particle density psiT) for a potential strength of cr = 20 and several values of 
7 > 0. The shape of Ps{t) with finite 7 is similar to the shape of the curve with 7 = 0, see Q, but all surviving 
particle densities for finite 7 are smaller than for 7 = 0. The results of figure ([T]) confirm, that the only influence 
of the particle-interaction on quantum reflection is indeed given by a slight acceleration of the decay of the surviving 
particle density inside the trap, that stems from an enhanced total energy. Other effects due to particle-interaction 
should modify the shape of the densities as functions of time, but the typical plateaus, the genuine pattern of quantum 
reflection already found in [3[ are conserved. The values of table jlll show that the surviving particle fractions after 
the scaled time r = 1 - approximatly one half of a second for the alkali atoms, see Eq. Q - depend only weakly 
on the self-interaction when 7 < 1 . Especially for higher values of a, the surviving particle densities for r = 1 
and 7 < 1 give still good results for trapping. Compared with a freely spreading wave-packet, which has a value 
Ps(t = 1) = 0.005, 0, the surviving particle densities for 7 = 1 are still enhanced by factors 10 to 40. For 
7 > 1, the particle density inside the trap starts to decay rapidly. For this regime, the mechanism of quantum 
reflection, although still working, is not strong enough for an effective trapping of atomic wave-packets, because the 
total energy of the system is too strongly increased by the interaction potential. As seen from Fig. ([2]), the strong 
repulsive interaction dominates the system at the beginning of the time-evolution, but as the atomic wave-packet 
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Ps{t = 1), 7 = 


0.1 0.5 1.0 5.0 


a = 20 


0.11 0.10 0.071 0.056 0.019 


(7 = 30 


0.19 0.17 0.14 0.11 0.043 


(J = 40 


0.27 0.24 0.20 0.16 0.06 


cr = 50 


0.33 0.31 0.25 0.19 0.09 



TABLE II: Listed are the surviving particle densities for several values of a and 7 for repulsive interaction. 

a = 20, Y=5 
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FIG. 2; Shown are the total energy of the system £'tot(T) (full line), the kinetic energy fkinCT") (dashed line) and the repulsive particle- 
interaction Eii^tir) (dotted line) for cr = 40 and 7 = 5. 

evolves in time, there is a continous current density of high-energy components beyond the edge of the step and a 
change of the shape of the wave-packet due to its motion on the step. The change of the shape influences the self- 
interacting potential and continously reduces its influence on the dynamics. But as the interaction energy decreases, 
it is transformed into kinetic energy. This leads to a faster motion and thus to a reduction of the effect of quantum 
reflection when large fractions of high-energetic components reach the edge of the step. Already after a time t « 0.1 
the system has strongly cooled down, but along with a strong loss of particle density. In short: the larger the kinetic 
energy of the wave-packet, the weaker is the effect of the quantum reflection. 

Therefore, we may conclude that a scaled particle-interaction strengths 7 < 1 effectively traps atoms. This is 
particulary true for large values of cr, which requires a weak potential strength /34, see Eqs. 

B. Dynamics with attractive particle-interaction 

In cases where the coupling constant 7 is negative, the particle-interaction is attractive. The additional, negative 
energy due to the self-interaction reduces the total energy of the system. As a consequence, the decay of the surviving 
particle fraction inside the trap is remarkably decelerated. The presence of attractive particle-interaction leads to a 
self-trapping of all lower momentum components which therefore are much more unlikely to even reach the edge of the 
trapping potential. In figure ([3]) we show surviving particle densities for different strength 7 of the self-interaction. 
For small interaction strength, 7 = —0.1, the dynamics is similar to the non-interacting case, 7 = 6, while for 
increasing interaction towards the critical interaction strength for collapse, we observe that the plateau structure 
becomes significantly washed out, 7 = —0.5, 7 = —0.62. 

Numerically we have determined the critical value 7c = —0.627. For attractive interactions stronger than 7c, the 
wave-packet dynamics will undergo a collapse. As will be discussed below, the value of 7c is universal for of the 
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FIG. 3: Surviving partile densities for tr = 40 and values of 7 = 0, —0.1, —0.5, —0.62 from bottom to top. The presence of the attractive 
particle-interaction increases the trapping effect of the quantum reflection. 

quantum reflection trap and a wider class of systems, which have no classical longitudinal (radial) confinement. 

The drastic change of the surviving particle density in the vincity of jc is reminiscent of the well-known phenomenon 
of a critical slowing down near the point of a phase-transition. The values of table pil[) show that the surviving particle 



Ps{t = 1), 7 = 


-0.1 -0.5 -0.6 -0.62 


cr = 20 


0.11 0.13 0.17 0.20 0.21 


cr = 30 


0.19 0.21 0.30 0.34 0.35 


0- = 40 


0.27 0.28 0.40 0.44 0.45 


cr = 50 


0.33 0.35 0.47 0.51 0.52 



TABLE III: Listed are the surviving particle densities for several values of a and 7 for attractive interaction. The reference time r = 1 
corresponds to around half a second for the here considered alkali species, see eq. Q. The surviving particle densities show that trapping 
by quantum reflection gives truely promising results in this regime. 

fractions psiT — 1) can achieve a strong enhancement when the particle- interaction strength 7 lies in the regime of 
> 7 > —0.627. Compared to the case of 7 = 0, the enhancement factor ranges from 1.6 to 2.0, but compared to 
the freely decaying wave-packet, ps{t = 1) = 0.005, the trapping mechanism achieves enhancement factors from 
42 for a = 20 to 104 for a — 50. Clearly, negative self-interaction, compared to the noninteracting case, significantly 
enhances the efficiency of the trap. Figure ^ illustrates the effect of trapping under self-interaction on the energies 
of the system as functions of time. At the beginning of the time-evolution, the total energy of the system is negative, 
and the particle-interaction dominates the system. The process of stabilization of the wave-packet by the self-trapping 
effect establishes a state of dynamical equilibrium between the kinetic energy and the interaction energy. When the 
total energy of the system has become positive, the wave-packet is no more in danger to collapse and has entered the 
regime were the kinetic energy, together with the quantum reflection, govern the behaviour of the system. However, 
the kinetic energy of the wave-packet has been reduced by the stabilization process, rendering quantum reflection 
even more effective. Together with figure ([3]) the behaviour of the energy clearly explains the trapping in this regime. 

C. Dynamical collapse for strong attractive interaction 

The dynamics of a wave-packet in the quantum reflection trap becomes unstable against collapse, when the in- 
teraction strength 7 becomes smaller than the critical value jc = —0.627. As the wave-packet evolves in time, the 
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a = 40, Y = -0.62 




FIG. 4: We have drawn the total energy of the system £tot{T) (full line), the kinetic energy £kin(''") (dashed line) and the repulsive 
particle-interaction £int(''") (dotted line) for cr = 40 and 7 = —0.62. It can clearly be seen how the energetic components of the system 
balance each other. The dynamical equilibrium between the kinetic energy and the interaction-energy stabilizes the system. 



influence of the attractive particle- interaction grows strongly, leading to a strong localization of the wave-packet. 
Strong localization goes along with high momenta, which lead to the destruction of the coherence of the wave-packet. 
When the wave-packet undergoes a complete disruption, it collapses. See the full curve in figure ([5]). The collapse 
goes along with a strong loss of particle density, that takes place during an extremely short interval of time. After the 
collapse of the wave-packet the attractive self-interaction is also destroyed and the effect of quantum reflection again 
dominates the behaviour of the system. The effect of quantum reflection manifests itself by stopping the evaporation 
and forcing the wave-packet to reshape. Thus, the remaining but strongly reduced density in the trap decays regularly 
as the system continues to evolve in time. Figure ^ shows the behaviour of both energies. To help the readers' eyes 
in separating the critical peaks of both energies, we have drawn two vertical lines. It can be clearly seen how the 
negative interaction energy grows to strong negative values during the time-evolution of the system. The kinetic en- 
ergy, on the positive half of the ordinate also grows during collapse in the attempt to achieve a dynamical equilibrium 
that stabilizes the system. However, by this peak in the kinetic energy the mechanism of quantum reflection becomes 
meaningless and a large fraction of particle density simply evaporates. Even after its peak the kinetic energy of the 
system is large enough to continuously drive a considerable fraction of particle density to the outside, such that the 
system experiences a self-cooling and the kinetic energy falls off again. When this has happened the effect of quantum 
reflection regains the control over the dynamics. Note, that the particle-interaction, due to the strong loss of particle 
density, remains negligible. 



IV. VARIATIONAL STABILITY ANALYSIS 



The stability of our quantum reflection trap can be accessed by variational techniques. A very elaborated and 
general variational approach to cold atom systems was given in [20|, where the critical interaction-strength jc for 
an anisotropic harmonically trapped 3D BEC-system in the absence of a longitudinal confinement was determined 
to 7c = —0.6204. The same problem is considered in [2l| numerically, where 7c = —0.627 was found. For our 
system, we have found the same critical value 7c = —0.627, despite a longitudinal (radial) confinement due to the 
atom-surface potential. The difference can be explained by the fact that the atom-surface potential is confining with 
respect to quantum reflection, whereas a harmonic oscillator potential is conflning due to classical reflection at the 
potential surface. When a longitudinal confinement is present in a harmonic oscillator system, the critical value is 
reduced to 7'^° = —0.57, [l^]. The absence of a (classical) longitudinal confinement due to a potential surface, along 
with the presence of a quantum confinement provided by atom-surface potentials thus enables the storage of a higher 
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FIG. 5: Surviving partile densities for cr = 40 and values of 7 = —0.63 (full line), 7c = —0.627 (dashed line) and 7 = (dot-dashed 
line) for comparison. The collapse of the wave-packet is indicated by the sudden fall-off of the surviving particle density. 




FIG. 6: Shown are the kinetic energy £'kin(''") (full line) and the repulsive particle-interaction £'int(T) (dashed line) for cr = 40 and 
7 = —0.63. The peak of the kinetic energy is responsible for the break-down of the trapping mechanism and the strong and sudden loss 
of particle density, see Fig. JSJ. 



number of atoms for attractive particle-interaction on macroscopic time-scales. 

To analyze the stability of our system, we will refer to a much simpler variational technique than the one suggested 
in [20|, see e.g. [1^, where the authors restrict themselves to a simple gaussian trial-function. Our choice for a radial 
trial-function defined on the whole space is 



N X exp 



X 



(11) 
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where N is the normahzation and a is the width. With Eq. (|lip we obtain a parametrized energy 



nia,a,j) = (^cj^^H^cj,) ^ ^ - exp [^a-^] - a eifc [a-'] j + -f^^ . (12) 

With Eq. we can estabihsh a relation between the particle-interaction 7 and the width of the gaussian state a 
by demanding that the variation of Ti. with respect to a vanishes, 

Snia,a,j{a)) = . (13) 

In Eq. (fT3|) we have kept cr as a parameter. Solving Eq. (I13|) for 7(a, cr) leads to 



7(a,cr) = - a + 



3 

2 a cr^ exp [-a^^] + ct^ (l - erf [a^^] ) - 



(14) 



Equation (jl4p defines a running coupling constant as a function of the initial width a of the wave-packet and the 
strength of the atom-surface potential a. The coupling constant depends implicitly on the initial kinetic energy of 
the wave-packet, which, as we have seen above, is crucial for the efficiency of the trapping mechanism. With a — 5 
in Eq.® we obtain a w 0.31 for Eq. (fTT|) . 

From Eq. (|14p. we obtain the phase-diagram of the quantum reflection trap. From the phase-diagram. Fig. ([7]), 
which we have drawn for all values of a considered above, the regions of stability can be read off easily and the value 
7c = —0.627 emerges as universal property of the quantum reflection trap. The horizontal and vertical lines mark 
our numerically determined value 7c = —0.627 and our initial data a = 5, respectively. 

As can be read off from figure ([7]), the simple gaussian trial- function delivers a critical region close to the exact 
value, such that we expect only small corrections from the general method of 20]. Also, our initial choice for a — 5, 
which matches the experimental conditions due to Jj, i2|] is justified as optimal choice. The universal region of the 
system is located left from the vertical line denoting a = 5 a — 0.31, where no diversification of 7 due to the 
potential strength a occurs. By increasing a, which means reducing a, the storage capacity of the system is reduced, 
because the critical value 7c is not accessible anymore. Vice versa the same occurs, because by decreasing a, which 
means increasing a, the region of universality is left and the stability of the system is only granted along the phase- 
lines described by the strength of the atom-surface potential tr, requiring a higher value for 7, which also reduces the 
storage capacity of the system for attractive particle-interaction. 

Lastly, Fig. ([7]), together with our numerical results from above proves, that the system shows not much sensitivity 
to the shape of the initial state, which neccessarily includes some arbitrariness in theoretical considerartions. 



V. SUMMARY AND CONCLUSION 



We have investigated the possibility of trapping cold atoms solely by the mechanism of quantum reflection. Our 
investigations have been carried out with the simplest possible model of a quantum reflection trap, a sphere. Typically 
the radius L of such a trapping sphere is of the order 10^ [a.u.]. The quantum reflection is mainly controlled by the 
strength of the atom-surface potential a = L/fi^. The larger cr, the better does quantum reflection work. The 
strength parameter a can be controlled either by the spatial extension of the trap, L, or by the atomic interaction 
parameter /34. Small values of (34 can be achieved, e.g., by using dielectric instead of perfectly conducting surfaces, 
see [ill [1^ and references therein. However, instead of controlling cr by surface-engineering, it may be much easier to 
control the value of 7 by a Feshbach-resonance. 

The inclusion of repulsive particle interaction, 7 > depletes the surviving particle density inside the trap in 
comparison to a system evolving without particle interaction. As long as the interaction strength 7 is smaller than 
unity, the depletion of the surviving particle density due to the repulsive interaction can be compensated by increasing 
the atom-surface interaction strength cr, see table (HI)) . 

The best results for trapping atoms by quantum reflection are achieved, when the particle interaction is attractive. 
The system remains stable when 7 lies in the range of > 7 > 7c. The critical value 7c — —0.627 is a universal 
property of our quantum reflection trap. This value of 7c makes clear that cold atom systems confined by an 
atom-surface potential belong to a wider class of systems with classical longitudinal freedom. Classical means, that 
no confinement due to a potential surface is present. The mechanism of quantum reflection acting in the case of 
atom-surface potentials establishes a quantum confinement. The quantum confinement has the advantage that more 
particles can be stored as in the case of a classical confinement, where the critical value is 7c° = —0.57. 

Approaching the critical value 7c from above, the typical plateau structure of ps{t) is washed out into slowly varying 
density oscillations, see Fig. ([3]). Close to the criticality, this behaviour is reminiscent of the well-known phenomenon 
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FIG. 7; Phase-diagram of the quantum reflection trap obtained by a simple gaussian variational approach. The vertical line a = 0.31 <-» 
a = 5 marks our initial data, the horizontal line jc = —0.627 marks our numerical result for the critical interaction-strength. The region 
encircled by the ellipse is the critical region of the system where transitions between universal and non-universal behaviour, as well as 
transitions between stable and unstable states of the condensate occur. The region below cr = supports only unstable states. The arrows 
drawn on the phase-line of cr = 10 indicate exemplarily for any other value of tr, that a stable state for a certain value of a becomes 
unstable when its parametric set (a, 7) is crossing the corresponding phase-line from left to right and vice versa. 



of the critical slowing down near the point of a phase-transition. Critical slowing down goes along with the existence 
of long-ranged fluctuations of the systems internal modes. It may be these long-ranged fluctuations that stabilize the 
wave-packet against the influence of the attractive self- interaction above 7c. 

For attractive particle-interactions 7 < —0.627, the internal motion of the wave-packet is not capable to stabilize 
the sytem. As the system evolves in time the wave-packet suffers a collapse. The collapse heats up the system because 
the interaction-energy is almost completely transformed into kinetic energy; along with the heating up, the system 
experiences a sudden loss of large fractions of particle density. 

Our results clearly show, that the mechanism of quantum reflection remains a promising tool to trap cold atoms 
when particle interaction is included and tuned. However, as our variational analysis has revealed, the parameters of 
the system cannot be chosen arbitrarily. For best results, they must allow the system to evlove in the universal region 
close to 7c. 

For alkali atoms, there are realistic surviving particle densities up to 50 % for times around half a second. From 
table (HI, the most promising candidate for trapping should be lithium, where the particle-interaction strength is 
already attractive. 
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